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Transmission eigenchannels and associated eigenvalues, that give a full account of wave propaga¬ 
tion in random media, have recently emerged as a major theme in theoretical and applied optics. 

Here we demonstrate, both analytically and numerically, that in quasi one-dimensional (ID) diffusive 
samples, their behavior is governed mostly by the asymmetry in the reflections of the sample edges 
rather than by the absolute values of the reflection coefficients themselves. We show that there exists 
a threshold value of the asymmetry parameter, below which high transmission eigenchannels exist, 
giving rise to a singularity in the distribution of the transmission eigenvalues, p(T —>■ 1) ~ (1 —T)”®. 

At the threshold, p(T) exhibits critical statistics with a distinct singularity ~ (1 — T)“3; above 
it the high transmission eigenchannels disappear and p(T) vanishes for T exceeding a maximal 
transmission eigenvalue. We show that such statistical behavior of the transmission eigenvalues 
can be explained in terms of effective cavities (resonators), analogous to those in which the states 
are trapped in ID strong Anderson localization. In particular, the p(T)-transition can be mapped 
onto the shuffling of the resonator with perfect transmittance from the sample center to the edge 
with stronger reflection. We also find a similar transition in the distribution of resonant transmit- 
tances in ID layered samples. These results reveal a physical connection between high transmission 
eigenchannels in diffusive systems and ID strong Anderson localization. They open up a fresh op¬ 
portunity for practically useful application: controlling the transparency of opaque media by tuning 
their coupling to the environment. 

PACS numbers: 42.25.Dd, 42.25.Bs, 71.23.An 


I. INTRODUCTION 


Wave progresses in random media exhibits rich 
physic^iM^, which finds numerous practical applications 
ranging from electron devices to optical communications 
and imaging. When a wave is incident on an open ran¬ 
dom medium it is decomposed into a number of “partial 
waves” that propagate independently along natural chan¬ 
nels - the so-called transmission eigenchannels - and are 
superposed again when they leave the medium. These 
channels can be obtained from the transmission matrix 
t, whose elements are coefhcients of field transmission 
through random media. The singular value decomposi¬ 
tion of this matrix, t = gives the wave¬ 

forms at the input and output edges of the nth trans¬ 
mission eigenchannel, i.e., the unit vectors u„ and v„, 
respectively and the corresponding transmission eigen¬ 
values rJS. The transmission eigenvalues and the eigen¬ 
channel structure - the corresponding spatial profiles of 
the wave fields inside the medium - give a full account 
of wave propagation in the interior of open random me¬ 
dia. In recent years, the coherent control of the inci¬ 
dent classical wave field has made it possible to control 
transmitted waves (e.g., Refs. EHUl), which has substan¬ 
tially advanced studies of wave propagation in random 
media. Subsequent investigations have been extended 
from the traditional subject of global transport behavior 
(e.g., transmission from the input to output edge) to the 
new realm of the eigenchannel structure (see Ref. |T3 for 
a review). 


Edge reflection arises from the refractive index mis¬ 
match at the surface of a sample, and is ubiquitous in all 
dielectric materials. The importance of edge reflection to 
the study diffusive wave transport was first pointed out 
in Ref. [TSl However, the investigations of its impact have 
so far largely been restricted to the average transmission 
and intensity of waves (see, e.g.. Refs. [T8U2T] as well as 
Ref. [5^ for a review). On the other hand, by reinjecting 
radiation that arrives at the edge from the interior of the 
sample, edge reflection strongly affects the coherence of 
waves and influences dramatically the transmission eigen¬ 
channels and eigenvalues at each individual realization. 
This fundamental issue yet remains largely unexplored. 

In a recent studji^ of wave transport through quasi 
one-dimensional (ID) diffusive media, it was discovered 
that an interesting “phase transition” occurs as the re¬ 
flection of one sample edge increases while the other edge 
remains transparent. This phase transition is seen in the 
distribution of transmission eigenvalues (DTE), defined 
as p(T) = ~ '^n)) with (•) being the disorder 

average. In the absence of edge reflection, 

where L is the sample lengthP^and ^(^ L) the localiza¬ 
tion length. This bimodal distribution was obtained in 
eighties of the past century, by Dorokho\E^ and Mello, 
Pereyra, and Kumail^ (see Ref. nzi for a review on the 
early status of this distributi on) and has r ecently received 
considerable renewed interestP^^^ElIllIll! Jt lays a foun¬ 
dation for mesocopic physics of electrons and photons. 
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TABLE I: Main characteristics of various phases. 


phase 

maximum transmission eigenvalue 

asymptotic behavior of p(T) 

O 

Tmax — 1 

p(r^i)~(i-r)-^ 

c 

Tlnax < 1 

P{T > Tkiax) = 0 

Critical 

Tmax — 1 

p(r^i)~(i-r)-3 


The singularity of po{T) for T —t 1 reflects the pres¬ 
ence of high transmission eigenchannels that dominate 
wave transport!^. An important question is: would these 
eigenchannels be blocked by edge reflection? In Ref. [231 
it was found that if edge reflection is below a certain 
critical value, they still exist, and the DTE displays the 
same singularity as po(T). Above the critical value, per¬ 
fectly transmitting eigenchannels disappear and the DTE 
is unimodal. At the critical value, the system undergoes a 
sharp transition and critical statistics emerges, display¬ 
ing a distinct singularity. This exact result completely 
washes out the common belief, i.e., that the main ef¬ 
fect of edge reflection is to elongate the sample length 
as enforced ^the well-known result of the average total 
transmissiorP^"2S! obtained by u sing eith er the diffusion 
model or radiative transfer theory^SMEil Instead, it sug¬ 
gests that even for diffusive waves, edge reflection affects 
significantly properties of transmission eigenchannels and 
eigenvalues. 

Many fundamental questions are thereby opened up. 
(i) How universal is the DTE transition? In fact, the 
finding of Ref. [23] relies on an exact solution, obtained 
owing largely to the system’s simple (and somewhat ar¬ 
tificial) construction. That is, one sample edge is trans¬ 
parent and the other reflective. A question naturally 
arises: what happens to realistic systems where both 
sample edges are generally semitransparent? (Studies 
of the complex eigenvalues of the non-Hermitian Hamil¬ 
tonian corresponding to these systems have shown that 
asymmetry of edge coupling to the environment leads to 
interesting transport phenomeneP^.) Eor such systems, 
would the phase transition still occur and the universal¬ 
ity (e.g., the critical statistics of transmission eigenval¬ 
ues) be affected? (ii) The origin of the transition has so 
far remained unclear. Notably, whether and how does 
wave interference gives rise to this transition? 

The purpose of this work is to extensively explore 
these subjects. To this end we consider wave transport 
through quasi ID diffusive samples with arbitrary edge 
reflectivities. We show analytically and numerically that 
the asymmetry of the surface interaction at two edges 
of the sample leads to rich phase transition phenomena. 
Specifically, for the asymmetry parameter below a cer¬ 
tain threshold value - even though the reflections of both 
edges could be strong - perfectly transmitting eigenchan¬ 
nels are opened, and the DTE exhibits a (1 —T)”® singu¬ 
larity (dubbed “0-phase”). Above the threshold, these 
channels are closed (dubbed “C-phase”) and the DTE 
becomes unimodal, vanishing above the maximum eigen¬ 
value Tjnax (< !)■ At the threshold, critical statistics 


with a (1 — T)“3 singularity occurs (dubbed “Critical 
phase”). The main properties of these three phases are 
summarized in Table [T] As we will show below, these 
asymptotic behaviors are insensitive to the system’s de¬ 
tails and thereby universal. Our analytical and numeri¬ 
cal analysis suggest that the perfectly transmitting (with 
the eigenvalue T = 1) eigenchannel is associated with 
an effective resonator bound to the center of the corre¬ 
sponding eigenchannel profile. Quite surprisingly, these 
resonators, although exist in diffusive samples, are of the 
same physical nature as the effective cavities, in which 
the eigenstates of the strongly disordered ID systems are 
localizecP^. Much as in the case of strong localization, 
the effective resonators in diffusive samples are formed 
by disorder-induced barriers and are of a small size - of 
order of the transport mean free path. When the asym¬ 
metry parameter increases, the centers of the eigenchan¬ 
nel profile and the resonator move towards the sample 
edge with stronger reflection. As they arrive at the edge, 
the perfectly transmitting eigenchannel disappears and 
the DTE transition occurs. 

We remark that the bimodal distribution 0 has been 
derived by vari ous theoretical methods in the past three 
decadeP^^^^^^lHall very recently received experimen¬ 
tal conflrmatiorfi^. In this work, we show that, physi¬ 
cally, the square root singularity in this distribution can 
be attributed to that the resonators are homogeneously 
distributed inside the sample. 

The remainder of this paper is organized as follows. In 
Sec. |nj we give a brief digest of the main results, and 


qualitatively discuss their physical meaning. In Sec. HI 


we present an analytical theory for DTE of quasi ID dis¬ 
ordered systems. In Sec. |IV[ a numerical verification of 
the analytical results is given. In Sec.jVj we introduce the 
ID resonators model, and use it to explain the physics 
of statistical behavior of the transmission eigenvalues in 
diffusive samples of higher dimension. We conclude and 
briefly discuss the results in Sec. in Sec. Ml Some tech¬ 
nical details are given in Appendices [^|D 


II. SUMMARY OF MAIN RESULTS AND 
THEIR PHYSICAL MEANINGS 

A. Phase structure of DTE 

We note that compared to the quasi ID sample with¬ 
out edge reflections, two new length scales appear in the 
presence of edge reflections. That is, the characteristic 
length Zb (z(,) associates with the reflection of the right 
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FIG. 1: Main panel: the phase diagram of DTE. Red lines 
correspond to the critical phase (|C ~ C^l = !)• Dashed line 
I (II) demonstrates the single (double) DTE transition when 
one edge reflection (i.e., ^') is fixed and the second (i.e., () 
increases. Inset: Macroscopic interpretation of phase param¬ 
eters C, The total resistance (inverse average conductance) 
of the sample is the sum of the left and right edge resistance 
and the bulk resistance. They are , ( and 1, respectively 
when rescaled by the bulk resistance L/^. 


(left) sample edge. More precisely, the ratio of this length 
to the transport mean free path, £, depends only on the 
reflection of the corresponding sample edge and mono- 
tonically increases with this edge reflection: it is of order 
unity when edge reflection vanishes and diverges when 
edge reflection reaches unity. The explicit expression of 
this ratio has been obtained by various methods over sev¬ 
eral decade^^ and, most recently, by the field theoretic 
approactP^. Because it is not essential to the present 
work, we shall not discuss this further. 

By using first-principles field theory, we find that the 
DTE or, equivalently, the factor 


/(r) = 


p{r) 

PoiT) 


( 2 ) 


characterizing the deviation from the bimodal distribu¬ 
tion po(T), depends on only two dimensionless parame¬ 
ters, C = ZhjL, = z'^jL. The physical meaning of these 
two parameters will become clearer in Sec. IIB Here 
we emphasize that the details of sample structure (e.g., 
disorder configuration, edge reflection, etc.) only affect 
their values as well as the average bulk conductance ^/L, 
namely, the conductance of the sample without the end 
reflections. The parameter of is irrelevant since it is 
an overall factor of po{T). So the behavior of DTE is uni¬ 
versal. For 1C —< 1, the high transmission eigenchan- 
nels are opened giving rise to a singularity ~ (1 — T )“2 


as T —>-1. For |C — Cl = 1, although the high transmis¬ 
sion eigenchannels are still present, a distinct singularity 
~ (1 — T)“s occurs. For |C — Cl > !> ttie high trans¬ 
mission eigenchannels are blocked; correspondingly, the 
asymptotic behavior of DTE undergoes a drastic change: 
p{T —>■ 1) = 0 and consequently p{T) becomes unimodal 
(with the peak near zero). These asymptotic behaviors 
as T —t 1 do not depend on the specific values of C,C^ 
rather on whether the asymmetry parameter |C — C^l is 
smaller (larger) than or equal to unity. This feature al¬ 
lows one to define three phases - O, C, and Critical - 
where the DTE behaves in qualitatively different ways 
(see Table |T] for a summary), resulting in the phase 
diagram (Fig. [^. It is symmetric with respect to the 
line of C = C^ reflecting the invariance of the DTE with 
respect to the exchange of the two edge reflections. 

Whenever a path in this phase diagram crosses a criti¬ 
cal line, the high transmission eigenchannels are switched 
on or off. As such, tuning the values of C,' leads to rich 
phase transition phenomena. To show this, we consider 
a path corresponding to fixing the reflection of one edge 
(say, the left, i.e., (') while increasing that at the other 
(i.e., C). We find that the DTE exhibits either single 
or double transitions: single transition occurs for < 1 
(Line I in Fig. [^; double transitions occur for C' > 1 
(Line II). We see that the previous result!^ corresponds to 
the special line of C' = 0 (or ^ = 0) in this phase diagram. 
Importantly, the stripe-like regime corresponding to the 
O-phase extends to infinity. This implies that the high 
transmission eigenchannels can still exist even when edge 
reflections are large, provided that they are not strongly 
asymmetric. This phenomenon resembles resonant trans¬ 
mission and signals strong interference origin of the high 
transmission eigenchannel. We shall explore this below. 


B. Physical meaning of parameters 

To better understand the physical meaning of con¬ 
trolled parameters we consider the average conduc¬ 
tance g defined as 


g= f dTTp{T). (3) 

Jo 

From the analytical theory developed below for p(T) it 
follows that, for arbitrary strength of edge reflections. 


L -\- zi) z'f^ 


( 4 ) 


We are not aware of any exact derivations of this result 
in the literature for arbitrary values of z^jL and z(,/L, 
although Eq. Q has been used in Ref. [TTJ Equation (|^ 
gives Ohm’s law (inset of Fig. [^: the reflection of the left 
(right) edge introduces an edge resistance (-^fc/C) in 
series with the bulk resistance (namely the resistance of 
the sample without edge reflections) L/^. This shows 
that, in sharp contrast to DTE, g does not exhibit any 








4 


criticality. Moreover, it shows that (, (' are the ratios of 
the corresponding edge resistance to the bulk resistance. 

Therefore, the role of these two parameters is two-fold: 
at the macroscopic level (i.e., as macroscopic observables 
such as the average conductance are concerned), they 
are essentially the edge resistors in series with the bulk 
one; at the mesoscopic level (i.e., as mesoscopic quantities 
such as the DTE are concerned), they play the role of 
phase parameters. 


C. Physical picture 



The transmission eigenchannel is associated with an 
energy profile (integrated over the transverse coordinate) 
across the sample. This profile is denoted as Wrix), 
with X being the distance to the left sample edge and 
T the corresponding eigenvalue. For perfectly transmit¬ 
ting eigenchannels (T = 1) in diffusive samples without 
edge reflections, this profile is a parabola with its cen¬ 
ter at X = Ll2 (Fig. red curve in the upper paneljP. 
In the presence of edge reflections, this profile remains 
parabolic but its center moves to the side with stronger 
edge reflection, i.e., larger edge resistance (Fig. blue 
curve in the upper panel). As we will see later, analyti¬ 
cal and numerical analyses suggest that, strikingly, there 
is a ID resonator bound to the profile Wq-^iix) (Fig. 
lower panel). This resonator is small, with a size of or¬ 
der of the transport mean free path, 1. It is formed by 
disordered barriers. It is analogous to an effective cav¬ 
ity in which the states are trapped in the regime of ID 
strong localizatioiP^. Importantly, it gives rise to perfect 
transmission (but not to the eigenchannel profile which is 
essentially the probability density for a diffusive wave to 
return to a cross section at depth x). If the asymmetry in 
edge reflections is increased continuously, the resonator 
moves together with the center of profile W-j-^ilx) to the 
sample edge, and eventually the DTE transition occurs. 


III. ANALYTICAL THEORY OF DTE 

In this section we study analytically impacts of edge 
reflection on the DTE. An essential difference from the 
earlier worlP^ is that in Ref. [52] the reflection appears 
only at one sample edge while in the present work both 
edges have finite reflectivities. As we show below, this 
difference introduces even more interesting physical phe¬ 
nomena. The generalization of the earlier analytical ap¬ 
proach to the present problem is, however, a highly non¬ 
trivial task. 


FIG. 2: Effective resonator related to a transmission eigen¬ 
channel with T = 1 is located in the center of the diffusive 
sample (left lower panel, red area) and is of the size of the 
transport mean free path £ (right lower panel). This resonator 
is bound to the center of the corresponding eigenchannel pro¬ 
file Wq-^iix) (upper panel, red curve). When the reflections 
of the sample edges are asymmetric, the center of the pro¬ 
file is shifted towards the edge with stronger reflection (upper 
panel, blue curve), and so is the resonator center (left lower 
panel, blue area). 


below. Our starting point is an exact expression for the 
DTE, 


P{T) 

n^) 


1 

27r 


[F{(j) -I- in) + F*{(j) + in)] 


d(j) 


1 

2 


sinhc^i 



ttf 

1 + sinh^((/)/2)ttt 


( 5 ) 

( 6 ) 


where (j) is understood as 4> — iS with S a positive in¬ 
finitesimal. The transmission eigenvalue F is related to 
the parameter (j) through 

T = cosh-^{(j)/2). (7) 

Then, it is a canonical method of casting the function, 
F((j)), into a functional integral over the supersymmetric 
field Q{x), where Q = {Q^a'} is a 4x 4 supermatrix, with 
A, A' = 1,2 denoting the advanced-retarded (‘ar’) sector 
and a, a' = f, b the fermionic-bosonic (‘fb’) sectop^S. Be¬ 
cause the results in this work do not depend on whether 
the time-reversal symmetry is present, we consider the 
system with broken time-reversal symmetry. This gives 
rise to a 4 X 4 supermatrix structure; otherwise one has 
to introduce additional matrix index to accommodate the 
time-reversal symmetry. The result reads 


A. General formalism 

We adopt the same general formalism as that of 
Ref. 1231 with an importance difference that we outline 


F(4>) = 


n(2zbQc»,,Q-[Q,r])U=i=0 


D[Q] 


(2z'QO,Q-S[Q.A])U=o=0 
^iQd.Q)liU=o,e=i^e-l ^0 (g) 
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with ‘str’ being the supertrace. Here A and F are con¬ 
stant supermatrices, 


A = 



( 9 ) 


r 


/ cos 9 
I i sin 6 


—i sin 0 \ 
— cos 6 J 


cosh(/) sinh0 Vin'i 
— sinh (j) — cosh cj) A ' 


with 0 < 9 < TT. Most importantly, edge reflection 
imposes boundary constraints on the supermatrix field, 
2z'f^QdxQ + [Q, A] = 0, at the left edge and 2zbQdxQ — 
[Q,r] = 0 at the right. As first shown in Ref. |40l these 
boundary constraints are the only effect of sample open¬ 
ness, namely the exchange of wave energy of media with 
external environments on the edges, on field theory con¬ 
structions. This notwithstanding, the constraints take 
rich physical wave phenomena in open media into ac¬ 
count (for a review, see Ref. HH). 

Equations ([^, and (|^ constitute an exact for¬ 
malism for calculating the DTE. We emphasize that this 
formalism is very general and valid for both quasi ID dif¬ 
fusive and localized systems, although the detailed treat¬ 
ments are very different. Also, the values of Z},, could 
be arbitrarily large: this advantage allows us to thor¬ 
oughly explore impacts of arbitrarily strong edge reflec¬ 
tions. In the remainder of this section we apply it to the 
former case where L 



FIG. 3: Deviation factor / as a function of the edge reflec¬ 
tion parameter ^ is fixed). DTE exhibits a single tran¬ 
sition C —>■ Critical —>■ O as increases from zero and 
is fix ed t o 0.6. Th e deviation factor is obtained by solving 
Eqs. (291 and (1301) numerically. The transition corresponds 


to (^ = 1.6 with /(T = 1) = 0 (dark blue line): above this 
value the probability density of perfectly transmitting eigen- 
channels is fully suppressed. From top to bottom the value of 
C is 0,0.3, 0.6, 0.9,1.2,1.6, 2, 3, and 5. 


B. Saddle point configurations 


Because of 1 the functional integral over the 

Q-field is dominated by fluctuations around the saddle 
point configurations which satisfy 

dxiQdxQ) = 0 ( 11 ) 


and the boundary conditions 


{2z'ijQdxQ + [Q, A])|2,^o — 0, 

{2zbQdxQ-[Q,mx=L=0. (12) 


The solution to the saddle point equation (11) has the 
same structure as T given by Eq. (10), i.e.. 


Q{x) 


f cos0(a;) 
jsin0(a;) 


—i sin 0(a;) 
— cos 0 (a:) 


ar 


f cosh$(x) 
— sinh<i)(a:) 


sinh $(a;) 
— cosh$(a;) 



(13) 


With the substitution of Eq. (HI, the saddle point equa¬ 
tion 0 is reduced to 

dl^ = die = 0, (14) 

and the boundary conditions ( |I^ to 

(z(,5a;$ - sinh$) |a;=o 
= (z(,5a;0 - sin0) | 3 ,=o = 0 (15) 


and 


{zbdx^ + sinh($ - (f))) \x=l 
= (^ 692 , 0 -hsm (0 - 0 )) = 0 , (16) 


respectively. 


The solution to Eq. (14) has the general form of 


<i)(a:) = C^x/L + (j)o, 
Q{x) = CexiL + 9o, 


(17) 

(18) 


where the coefficients G R, ((>0 > 0 and 0 < 9o < tt. 
Due to the compactness of the fermionic component, i.e., 
the 27r-periodicity of the sine and cosine functions, there 
are a family of saddle point solutions, each of which 
corresponds to a saddle point action proportional to 
|;(C'e^i 0 + C^). Assuming that the compactness does not 
play any role, we may carry out the analytic continua¬ 
tion of the second equations in the boundary constraints 
(15) and (16), which gives Ce^icj) = iC^. This is the sad¬ 
dle point configuration with the smallest action which is 
zero. It preserves the supersymmetry. The other sad¬ 
dle point configurations arise from the compactness and 
break the supersymmetry. They are gapped by an ac¬ 
tion ~ C>(|;)- This action is large for diffusive waves and 
the supersymmetry broken saddle points are can thereby 
be neglected. Furthermore, we may ignore fluctuations 
around the saddle point in the pre-exponential factor of 
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TABLE II: 


Numerical values of a and /3 obtained by numerically solving Eqs. (291 and (301 at = 1. 


C 1 

3 

5 

7 

9 

11 

13 

15 

a 0.333331 

0.333228 

0.333078 

0.332906 

0.332716 

0.332510 

0.332292 

0.332062 

P 0.333834 

0.339621 

0.342927 

0.345654 

0.348075 

0.350291 

0.352356 

0.354303 


Eq. ^ since they only give rise to corrections of lower 
order. Finally, by integrating out Gaussian fluctuations 
around the supersymmetric saddle point, which gives a 
functional superdeterminant of unity aA, 0 = i(p because 
of the supersymmetry, we reduce Eq. to 

F{4>) = -I^C4,. (19) 


This is uniquely determined by the coefficient and we 
therefore focus on the solution of <i)(a:) below. 

For the limiting case of Zh, z'^ — 0, the boundary con¬ 
straints (15) and ([T^ [for $(a;)] are reduced to 


$ 




= 0 , $1 


x—L — 


( 20 ) 


Combined with Eq. ([T7| ) they give $(a:) = 4>x/L. Substi¬ 
tuting it into Eqs. ([^ and (19) gives Eq. ([^. 

By using Eqs. (§ and ( |19[ ) and taking into account the 
parametrization of T we obtain 


/(<()) = - C^.„)l{2m) = AC^/{2i7r). (21) 

According to this equation we need in principle to find 
first and then make analytic continuations. But this is a 
difficult task because as we will see shortly, the equation 
satisfied by is transcendental. To overcome this dif¬ 
ficulty we will establish a closed equation of below 
and solve the equation. 


C. DTE transition 


From Eqs. (|15|) and (16) we find 


(C'G^ - sinh(/io) U=o = 0, 

(CG^ -b sinh(G0 + (jio - <(>)) \x=l = 0. 

Substitution Eq. ( [^ into Eq. ( [^ gives 

C' cosh(G0 - (j))C^ + (C 4 , 

cosh(()o =- 

smh(G,^ - 0) 

Combining Eqs. (22) and ( p4| ) we find 
[C cosHC^ - + CC^]^ 


- eel = 1. 


sinh^(G0 - 0) 

This equation may be rewritten as 

sinh^(G0 - 0) 


Ci = 


( 22 ) 

(23) 


(24) 


(25) 


2CC'cosh(G0-0) + C2 + C'"’ 

which implicitly gives Ctj, as a function of 0, (C, and 00 


(26) 
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EIG. 4: Deviation factor / as a function of the edge reflection 
parameter 0 (0' is fixed). DTE exhibits double transitions 
C — >■ Critical — >■ O — > Critical — > C when 0 increases from 
zero and 0' is fixed to 4. The deviation factor is obtained 
by solving Eqs. ( |29| l and ( |30| l numerically. The transitions 
correspond to 0 = 3 (upper panel, purple dashed line) and 
0 = 5 (lower panel, dark green dashed line), respectively. 
/(T = 1) is positive for 3 < 0 < 5 and otherwise zero. 


We introduce = (G0+i^-|-G0_i,r)/2. Then, by using 
Eq. (26) we obtain 


sinh^ 0_l_ sinh^ 0_ 


2G0 + ^^ 0/2 - acosh0+ + 6 + ^osh0_”+6 


(27) 







































7 


and 


2C^AC^ = 


sinh^ ')/>+ 


sinn ■!/) 


a cosh -ip^ +b a cosh ip- + b ’ 


(28) 


where 'ip± = C^±iTr — (<(' ± ^tt), a = 2 CC', and 6 = 
^2 _|_ ^/2 ^ Equations ( [^ and (28) are equivalent to (see 
Appendix [A] for the derivation) 


=i2 2 ^ 2 ^ —a(cosh 2 (C '0 — </>) + cos27r/ — 2) cosh(C 0 — (p) costt/ + 5(cosh2(C'^ — (p) cos27r/ — 1) 


2{C.-7r^n = 


a^(cosh2(C0 — <(>)+ cos27r/)/2 — 2a6cosh(C'0 — iji) costt/ + 6^ 


(29) 


47rC'0/ 


—a(cosh2(C^ — <(')+ cos 27r/ + 2) sinh(C'^ — (p) siuTr/ + 6sinh2(C'0 — (p) sin27r/ 
a^(cosh2(C'^ — <(') + cos27r/)/2 — 2a6cosh(C'0 — (p) costt/ + 6^ 


(30) 


Recall that is real and / > 0. Equations (29) and (30) 
constitute the closed equations for and /. Giving /, 
we may solve them numerically and find the deviation 
factor /. The representative results of / are shown in 
Figs. and 1^ 


As shown in Fig. if the reflection of the left edge, 
i.e., C' 1 is fixed and sufficiently small, upon increasing 
the reflection of the right edge, i.e., the DTE exhibits 
a single transition at certain critical value of / with a 
hallmark of /(T = 1 ) = 0 : below the critical value the 
maximum transmission eigenvalue is smaller than unity 
while above the critical value it is unity. This transition 
is similar to the previous resulfP^ where is set to zero. 
As shown in that work, in such limiting case the transi¬ 
tion coincides (but generally not) with a DTE transition 
predicted for a completely different setup, i.e., a normal 
metal with a single barrier placed inside the sampl^^. 


The behavior is even more interesting when the fixed 
C is large. (For this case we are not aware of any analogs 
of the results derived below in other wave systems.) As 
shown in Fig.[^ upon increasing C, the DTE exhibits dou¬ 
ble transitions at two critical values of C,: below (above) 
the lower (upper) critical value the maximum transmis¬ 
sion eigenvalue is smaller than unity while between these 
two critical values is unity; the critical points have a hall¬ 
mark of /(T = 1 ) = 0 . 


D. Exact criterion for DTE transition 


To find the condition under which this equation has a 
real solution of > 0, we rewrite Eq. (311 as 




F{C^) = 2C^{h — acosh.C^) — (cosh 2(70 — 1) = 0. (32) 


Performing the Taylor expansion for F{C^) gives 


F{C^)=2{b-a-l)cl-J2 


n=2 


2a 


(2(n-l))! (2n)! 


c!” 


(33) 


for I ( 70 1 < 00 . From this we see that provided 

6 -a- 1 > 0^ IC-C'I > 1, (34) 


F{C^) is a non-monotonic function of ( 70 : it first in¬ 
creases from zero, then decreases and eventually decays 
as ~ when is sufficiently large. In this case, 

Eq. (32) must have a positive root. If the inequality (34) 
is not satisfied, then F{C^) monotonically decreases from 
zero and Eq. (32) has only a trivial solution of (70 = 0 
and, as a result, /(T = 1) must not vanish. 


The inequality (34) defines the two regimes of (7-phase 
in Fig. [2 It implies that the perfectly transmitting 
eigenchannel is blocked only if edge reflections are highly 
asymmetric. Otherwise, |C — C'| < 1 which defines the 
regime of (7-phase in Fig.[T] The transition occurs at 


IC-C'I = 1. 


(35) 


Next, we derive the exact criterion for the transition. 
Because a hallmark of the latter is /(T = I) = 0, as 
mentioned above, we set (p = 0 [recalling that T = 
cosh“^((()/2)] and / = 0 in Eqs. (29) and (30). We find 
that the latter equation is satisfied automatically and the 
former is reduced to 



cosh 2(70 — 1 
b — a cosh (70 


(31) 


This gives the two critical lines in Fig. From this we 
see that for < 1 single transition occurs upon increas¬ 
ing / from zero, and the critical point is {(' + !)■ This 
phase structure is represented by Line I in Fig. In the 
particular case of C' = 0 this result is in agreement with 
that found in the previous worlP^. For > 1 two transi¬ 
tions occur upon increasing / from zero, and the critical 
points are (/' ± 1). This phase structure is represented 
by Line II in Fig. 


























E. Criticality of DTE transition 


-2 


1. Critical statistics 


From the results above we find that the DTE peak near 
T = 1 diverges as 

p(r^ 1 )^( 1 -r)-5, ic-c'l<i, (36) 


since /(T = 1) 7 ^ 0. In this part we study the asymptotic 
behavior of p(T —t 1) at the critical line. To this end we 
solve Eqs. (29) and (30) numerically. We find that for 
(/) —?► 0 corresponding to T —>■ 1 the solution has the 
general form as follows (cf. Fig.[^, 


C^ - 4> (X 4>P. 


(37) 


The exponents, a, /?, for different values of are given 
in Table [Tl| (Because both the values C ~ C = ±1 lead to 
the same results, without loss of generality we consider 
C — C = 1.) These numerical values give a « /3 « 1/3. 
In Appendix [B| we show that a stronger relation, i.e.. 


a = 13 = 1/3 


(38) 


exists. Then, combining Eqs. (37) and (38) and taking 
into account that ~ (1 — T)^'^or (/ —)■ 0 , we find 


p(r^ 1 )^( 1 -r)-5, ic-C'l = i- (39) 


We see that the divergence of the DTE peak around T = 
1 changes whenever the critical line is crossed, associated 
with the emergence or disappearance of high transmission 
eigenchannels. In the special case of vanishing C or , 
such a change in the power of the divergence from I /2 to 
1/3, has been found previouslji^. 



- 8.0 - 7.5 - 7.0 -6.5 - 6.0 -5.5 - 5.0 

ln(0) 

FIG. 5: The numerical solution of fitf)) for different values 
of ( and satisfying ^ — C' = 1 (solid lines). From top to 
bottom C = 1,3,5,7,9,11,13,15. To guide the eye we also 
plot a dashed line with a slope of 1/3. 


it is due to that the field theory, in a somehow auto¬ 
matical manner, captures collective modes such as the 
diffuson (see Ref. im for a pedagogical review). To see 
this, in Appendix [C] we investigate the Gaussian fluctua¬ 
tions in the field theory ^ and find the diffuson explic¬ 
itly. The result entails the canonical physical meaning of 
Zb, z'^ namely the so-called extrapolation lengtlP^. 


IV. NUMERICAL TEST OF DTE TRANSITION 


2. Average conductance 


Substituting p{T) = /(T)po(T) into the definition 
(§, we find that the average conductance, upon being 
rescaled by ^/L, is given by 


g X (L/o = 






(40) 


It depends only on the parameters (/ and C' ■ An interest¬ 
ing question is, does this observable exhibit any critical¬ 
ity? Because we cannot analytically calculate Eq. (40), 


we resort to carrying out the integration numerically. To 
this end we let (/, (/' vary over a wide range from 10 “^ to 
10^. As shown in Fig. the numerical results are well 
fitted by 1/(1 -I- C + CO- Therefore, the average conduc¬ 
tance is given by Eq. ([^, which is none but Ohm’s law 
and shows that the average conductance does not exhibit 
criticality. It agrees with the well-known resullP^^^ ob¬ 
tained byj^ng the diffusion model or radiative transfer 
theorySHMSi] This agreement is not accidental. Rather, 


In this section we put the analytic results derived in 
Sec. mil under numerical test. To this end we launch a 
scalar wave into a quasi ID disordered waveguide which 
is a 900 X 300 rectangular lattice and thereby locally two- 
dimensional. We simulate wave transport by using the 
standard recursive Green’s function methocP^HBl Specif¬ 
ically, we are interested in the Green’s function G(r, r') 
between grid points r = {0,y) and r = {L,y') at the left 
{x = 0) and right (x = L) edges, with y being the trans¬ 
verse coordinate. The lattice constant is the inverse wave 
number. The squared refractive index at each site fluc¬ 
tuates independently around the air background value, 
taking values randomly from the interval [0.7,1.3]. The 
wave velocity in the air background is set to unity. To 
create edge reflection, we add an additional layer of thick¬ 
ness 2 with uniform refractive index at the corresponding 
sample edge. For given values of refractive index at the 
left and right edges and disorder configuration, we nu¬ 
merically compute the transmission matrix t = {tba} in 
the basis of the empty waveguide modes, 4>a,b{y), where 
the indexes a, b are the labels of the left and right edges. 
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FIG. 6: Numerically calculated Eq. ( |40| | (dots and circles) 
shows that the average conductance obeys Ohm’s law irre¬ 
spective of the values of and . Note that the two data 
sets corresponding to Zb — 0.4, z), = 2 and Zb = 0.8, z'b = 1.6, 
respectively, are identical because they have the same value 
of Zb + z'b- The slope of the fitting solid red line is one. 


respectively. These matrix elements are given by 

ha = y/vbVa J dy j dy'(l)biy)'Pl{y')G{r,r'), (41) 

where Va is the group velocity of the empty waveguide 
mode a at the wave frequency. For simulations the chan¬ 
nel number is set to 100, i.e., t is a 100 x 100 matrix. By 
numerical diagonalization we find the transmission eigen¬ 
values, r„. Then, we repeat the same procedures for 6000 
disorder configurations. 


A. Ohm’s law 

With the transmission eigenvalues obtained by the 
numerical method described above, we compute the av¬ 
erage conductance given hy g = First of all, 

we set the refractive index of the dielectric layers at the 
left and right edges to unity so that no (edge) reflections 
arise. As shown in Fig. the resistance, g~^-, increases 
linearly with the sample length, L. This confirms that 
the sample is a diffusive conductor. Indeed, the slope of 
the straight line in Fig. is 6.33 x 10“®, and noticing 
g-^ = (L + 2zbo)/^ [cf. Eq. @] with Zbo being the ex¬ 
trapolation length in the absence of edge reflectior l^^F4 | 
we obtain the localization length ^ = 1.58 x 10"^ which is 
much larger than the sample length {^/L = 17.6). More¬ 
over, the intersection of the straight line with the vertical 
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500 600 700 800 900 

L 


FIG. 7: Numerical simulations show that in the absence 
of edge reflection the average conductance (squares) obeys 
Ohm’s law (bottom solid line). The slope of the line gives the 
inverse localization length, and its intersection with the 
axis L = 0 {2zbo/0 gives the value of Zbo- For different values 
of the refractive index. Ohm’s law (top and middle solid lines) 
remains valid, with the slope of the straight lines unchanged. 
The refractive indexes at the left and right edges are (1,1.9) 
(triangles) and (1,2.1) (circles), respectively. 


axis L = 0 gives 2zbo/^ = 9.52 x 10“^. With the substi¬ 
tution of the value of ^ we find Zbo = 75.31. 

For fixed different values of refractive indexes at the 
left and right edges, by varying the sample length and 
computing g we also confirm the general expression 0. 
Importantly, this result allows us to determine the value 
of the length Zb for different values of the dielectric con¬ 
stant e of edge dielectric layer. Specifically, we fix e of 
the dielectric layer at the left edge to be unity while 
vary that at the right. For each dielectric constant value 
corresponding to the right edge we compute the average 
conductance for different sample lengths. As mentioned 
above, g~^ increases linearly with L. For this straight 
line, we find that, as shown in Fig. the slope is inde¬ 
pendent of dielectric constant at the right edge, in agree¬ 
ment with the prediction of Eq. (|^, whereas the inter¬ 
section with the vertical axis L = 0 varies. Subtracting 
Zbo from the product of ^ and this intercept, we find Zb 
for corresponding e at the right edge. The results of Zb 
for different values of e are given in Table [Ill| 


B. DTE transition 


We now focus on the most interesting case of double 
phase transitions. To numerically explore the DTE tran¬ 
sition we fix the dielectric constant of the reflection layer 
at the left edge to 2.1 (corresponding to (^' = 4.60, see 
Table III) and tune that at the right edge from I to 2.4 


with the increment of 0.1. Because of > 1 the cor¬ 
responding phase structure must be the same as Line II 
in Fig. According to the exact criterion (35) and the 
numerical value of C given in Table III we predict that 
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TABLE III: The value of Zb obtained by numerical simulations. 


€ 

1.0 

1.5 

1.6 

1.7 

1.9 

2.0 

2.1 

2.2 

2.3 

2.4 

Zb{z'b) 

75.31 

92.01 

103.3 

162.5 

979.8 

2111 

4141 

7504 

12996 

21123 

C(C') 

0.08 

0.10 

0.11 

0.18 

1.09 

2.35 

4.60 

8.34 

14.44 

23.47 


the system undergoes double phase transitions when e in¬ 
creases and takes the values in this table: for e less than 
(or equal to) 2 . 0 , because of — C > 1 the system is in 
the C-phase (upper left regime in Fig. [^; for e = 2.1, 
because of — C = 0 the system is in the 0 -phase; for e 
larger than (or equal to) 2 . 2 , because of C — > 1 the 

system enters into the other O-phase regime (lower right 
regime in Fig. [^. These double phase transitions are ob¬ 
served numerically, as shown in Fig. where the colors 
of blue, red, and black correspond to two 0 -phases in 
different regimes of Fig. and the O-phase, respectively. 
(To find the critical phase is beyond our current numer¬ 
ical experimental reach since this would require the fine 
tuning of e at the sample edge.) 

Moreover, the analytic results of the /-factor (solid 
lines), with the values of given by Table III are 
in good agreement with simulation results without any 
fitting parameters. The smearing of the singularity is 
an effect of finite number of channels. In simulations, 
the transmission eigenvalue density at T is defined as 
the ratio of the number of eigenvalues in the interval of 
[T—|AT, T-I- 5 AT] (AT set to 0 . 01 ) to the total number 
of eigenvalues. Therefore, numerical experiments confirm 
that for sufficiently large e', so that /' > 1 , the system 
undergoes double DTE transitions as e at the right edge 
increases from unity. 

As we have analytically shown, the phase structure 
corresponding to Line I in Fig. where the single phase 
transition occurs, is similar to its limiting case corre¬ 
sponding to the line of /' = 0 in Fig. The latter has 
been confirmed in numerical simulations previouslj^^. So 
we do not discuss single phase transition here. 


V. PHYSICAL MECHANISM 

Above we have developed an exact analytical theory 
for the DTE transition and confirmed the transition by 
numerical simulations. In this section we further pro¬ 
vide a transparent physical picture based on combined 
analytical and numerical analyses. 



FIG. 8: Simulations confirm the occurrence of double phase 
transitions. The dielectric constant of the reflection layer at 
the left edge is fixed to e = 2.1 and that at right, e, increases. 
The simulation results are for e = 1.0 (blue empty circles), 1.9 
(blue solid circles), 2.1 (black solid diamonds), 2.3 (red empty 
squares) and 2.4 (red solid squares), respectively. The blue 
(red) circles correspond to the upper left (lower right) C-phase 
regime in Fig.[^where the high transmission eigenchannels are 
blocked; the black circles correspond to the O-phase where the 
high transmission eigenchannels are present. The simulation 
results are in good agreement with the analytical results of 
/(T) (solid lines) without any fitting parameters. 


also Ref. mi for a review). Based on this fact it is nat¬ 
ural to expect the eigenvalue, i.e., the parameter (j) in 
Eq. 0 depends implicitly on a certain position variable, 
X, associated with the eigenchannel structure. Loosely 
speaking, this variable characterizes how deep the en¬ 
ergy is pumped into the medium. Its physical meaning 
will become clearer below. 


1. (p-resonator center correspondence 


A. Mapping to resonator model 

We first note that, very recently, it has been 
discoverecP that the structure of eigenchannels, namely, 
the (spatial) profile of the wave energy density, is uni¬ 
versal. It is governed by the localization length asso¬ 
ciated with the corresponding transmission eigenvalu^^ 
and the (position-dependent) diffusion coefhcieniP^ (see 


Let us hypothesize (see Appendix |0 for more discus¬ 
sions) that the explicit dependence of (l){x) is 

(j) = 2{L-2x)/t (42) 

Recall that the left edge of the sample corresponds to 
a: = 0. Then, 

_ 9 L 2x ,..... 

T = cosh — j — (43) 
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• s=0.1 

■ s=0.5 

■ s=0.7 


with 
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Tr^t-'S** . . * 


FIG. 9: Simulations show that for sufficiently strong disorders 
(large s) DRT exhibits bimodal behavior. The samples are ID 
with transparent edges. 


and 


dT 


AL rVT^ 


= d w 


(!)^ 


(44) 


This shows that a bimodal distribution of T is equiv¬ 
alent to a uniform distribution of the position variable 
X. Indeed, the distribution given by the left-hand side 
of Eq. (441 differs from Eq. Q only in an unimportant 
overall factor. More precisely, it can be shown that pro¬ 
vided X is uniformly distributed over the sample, then 
the mapping: x ^ (j) leading to the bimodal distribution 
is unique (see Appendix]^ for the proof), and is given by 
Eq. (1^. 


According to Eqs. \A2\ and ( |43[ ), the eigenvalue T = 1 
corresponds to a: = L/2! In the presence of edge reflec¬ 
tion, Eq. (|^ suggests that effectively the sample length 
is extended by an amount of from the left edge and 
of zi, from the right. Taking this taken into account, we 
modify Eq. (431 as 


T = cosh 


-2 A- + Zb) — 2(x + z'l,) 
2 £ ’ 


and the eigenvalue T = 1 corresponds to 
x = {L + Zb - z'b)l2. 
Because 0 < a; < L, we find 

\zb - 41 < L. 


(45) 


(46) 


(47) 


This is identical to the inequality (34) for which perfectly 


transmitting eigenchannels are present. Therefore, the 


criterion for the DTE transition (35) is reproduced. 

To better understand the physical meaning of the po¬ 
sition variable x we note that the parametrization (43) 
has an equivalent expression. 


T{x) = 


AT'{x)T{x) 

[T'{x)+nx)f 


T\x) = e 




T(x) = 


(49) 


Surprisingly, Eq. (481 is identical to the expression for the 


transmittance through a disorder-induced ID resonator 
in ID random media in the localized regime, where the 
localization length is order of the transport mean free 
pathP^. [In fact, the expression is similar to the gen¬ 
eral expression ( |5^ below for resonant transmittance.] 
Because x,L — x are much larger than £, the resonator 
corresponding to T{x) given by Eq. (48) essentially is 
a point-like defect, placed at x and confined within a 
size ^ This analogy between the transmission 

eigenvalue and resonant transmittance suggests that the 
(high) transmission eigenchannels have a close relation 
to resonators and therefore are of strong interference ori¬ 
gin. The mapping of a high transmission eigenchannel 
in the system of higher dimension onto ID resonator 
makes sense perhaps because the eigenvalue T corre¬ 
sponds to the universal eigenchannel profile, and this 
profile is structureless in the transverse direction of the 
waveguid^. 


2. Analog of f-factor in resonator model 


As shown in Eq. (44), in samples without edge reflec¬ 


tions, the bimodal distribution follows from a uniform 
distribution of resonator centers. To find out what hap¬ 
pens in systems with reflective edges, we multiple both 
sides of Eq. (44) by the deviation factor / to obtain 


i f{T)dT 

2L rVT^ 


/[r(a:)]d(|). (50) 


The left-hand side gives the general expression of DTE, 
i.e., p{T) = /(T)po(T); the right-hand side shows that 
the physical meaning of the /-factor is the spatial den¬ 
sity of resonators. It then becomes clear that the closing 
of the perfectly transmitting eigenchannel can be inter¬ 
preted as the vanishing of resonators with T'{x) =T{x). 
Indeed, if the reflections of the left and right edges are so 
asymmetric that an effective resonator with T' {x) = T(x) 
cannot be introduced, then according to Eq. ( [4^ the 
maximum transmission eigenvalue must be smaller than 
unity. 


3. Normalization condition in resonator model 

We note that the bimodal distribution is not normal¬ 
izable. Because the integral po{T)dT suffers logarith¬ 
mic divergence due to the singularity for T —> 0. 
Therefore, we cut the integral at a certain exponentially 
small value, Tc « so that po{T)dT = N. On 


the other hand, from Eq. (48) it is easy to see that 


r(x) > Ttyp = 4e-"^/^ 


(48) 


(51) 
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FIG. 10: Simulations of ID samples show that the perfect res¬ 
onant transmittance (T = 1) disappears when the reflection 
coefficient (1 — Tmin) of one sample edge (the other is always 
transparent) exceeds a critical value. The disorder strength 
parameter s = 0.5. The vertical axis is the number of events 
and the distribution, therefore, is not normalized. 


FIG. 11: Simulations confirm that the critical value of the 
edge reflector r^rit, above which perfect resonant transmit- 
tances disappear, is the same as the maximal value Cmax of 
|rs(aj)| for the same sample without edge reflector. 


This implies that the cutoff of the bimodal distribution, 
Tc, has the meaning of the typical transmission, Tjyp, in 
the resonator model. 


B. Transition in distribution of the resonant 
transmittance (DRT) 

From analysis above we have seen that although the 
transmission eigenvalue and resonant transmittance are 
very different concepts, they have well-pronounced simi¬ 
larity. Naturally, we expect the DRT to exhibit a transi¬ 
tion similar to that of DTE. In this part we will address 
this issue. We perform numerical experiments on wave 
transmission through ID random media, namely layered 
samples in the strong localization regime. 


1. Numerical observations of DRT transition 

We randomly place a fixed amount (set to 50 in numer¬ 
ical experiments below) of scatterers in a ID chain, and 
let the reflection coefficients of scatterers be svi {i labels 
the scatterers). For each i, ri randomly takes a value 
from the interval of (—1,1). The constant s satisfying 
0 < s < 1 is independent of i. The disorder strength is 
controlled by the value of s. The dimensionless distances 
between nearest scatters are randomly distributed in the 
interval AS/Sq = 90%. We change the frequency 

w of incident waves in the narrow interval coq ± Aw, 
Aw/wo = 5%, and calculate the transmittance spectrum 
7~(uj) by using the standard transfer matrix approach. 
The numerical experiments are repeated for 10^ disor¬ 
der configurations. There are about 10^ resonant peaks 


(transmittance resonances) in the spectrum T(w) in any 
configurations in the given frequency interval, so that the 
total number of resonances analyzed below is about 10®. 

First, we study the case where both sample edges are 
transparent. We find that the background value of the 
spectrum T(w), namely the typical transmission coeffi¬ 
cient, is very close to zero in accord with the expression 
of Ttyp given in Eq. (51). The resonator corresponds to 
the local maximum in tnis transmittance spectrum, i.e., 
T = max{T(w)}. As shown in Fig. for sufficiently 
strong disorder strengths (so that the resonator centers 
are uniformly distributed in the sample) the DRT, de¬ 
noted as PriT), is in good agreement with the distri¬ 
bution ([^ while for moderate disorder strengths Pr[T) 
deviates from the behavior of Eq. 0 but is still bimodal. 

Next, we study effects of edge reflection on the DRT. 
For simplicity we consider only the case where one (say 
the left) edge is transparent while the reflection of the 
other edge, denoted as (1 — Tmin), increases. The sim¬ 
ulation results are shown in Fig. We see that when 
the reflection increases and exceeds certain critical value, 
the peak at T = 1 is fully suppressed, i.e., the perfect 
resonant transmittance (T = 1) disappears. When the 
internal reflection further increases, the highest resonant 
transmittance decreases. This behavior is the same as 
the single transition exhibited by the DTE for quasi ID 
disordered waveguides with one edge transparent and the 
other reflective. 


2. Physical mechanism for DRT transition 

We show below that the DRT transition can be well ex¬ 
plained by the resonant cavity theory for ID disordered 
system^. Suppose that the ID sample without edge re- 
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flection has a complex reflection coefficient rs(a;). This 
sample itself may be considered as a single semitrans¬ 
parent barrier, characterized by rs{uj)- Upon placing ad¬ 
ditional reflector after the output edge of the sample, 
which is described by a complex w-independent reflec¬ 
tion coefficient rend) we have two barriers that constitute 
a resonator. The total transmittance of this resonator is 


r(cd) 

Tsiuj) 


4rs(cd)Tniin 

(T,H+T„,in)'’ 

l-lr.Hp, T„,in = l-|rendr (53) 


This shows that the resonator can be perfectly transpar¬ 
ent only if |rs(a;)| = |rend| and the round trip phase shift 
is multiple of 2tt. When |rend| is small enough and the 
frequency varied in a rather broad range, there are many 
frequencies for which these conditions are satisfied. How 
does the number of such frequencies, W) evolve with 
increasing |rend I? Obviously, the number of frequency 
regions where |rs(a;)| = |rend| becomes smaller and 
decreases. When |rend| exceeds maximal value Cmax of 
|rs(a;)| (for given sample in given frequency range), the 
condition |rs(a;)| = |rend| cannot be satisfied by any fre¬ 
quency. Thus, Tmax IS the critical value of the edge re¬ 
flector above which perfect resonant transmittances dis¬ 
appear. 

To check this result numerically, we firs find Cmax for all 
ID random samples. Then, we add additional layer with 
reflection coefficient rend and select those T{uj) larger 
than 0.999 as perfect resonant transmittances. Then, 
we increase |rend| until no frequencies in the given range 
satisfy this selection rule. This value, denoted as rcrit) is 
the critical value for the DRT transition. As shown in 
Fig.[TTl simulations confirm rcrit = ''’max- 


C. Physical picture of DTE transition 

Based on the demonstrated similarities between the 
transmission eigenvalues and resonant transmittances, 
and between the DTE and DRT transitions, in what 
follows, we propose a physical explanation of the DTE 
transition. 

As shown in Ref. HI for quasi ID samples without edge 
reflections, the eigenchannel structure corresponding to 
perfect transmission is given by 

Wr=i{x) = l + Fiix), (54) 

where Fi{x) is essentially the probability density for a 
wave to return to a cross section at depth x in an open 
random medium. (We put the left edge as the input 
edge.) The independence of the eigenchannel profile 
W-r{x) on the transverse coordinate reflects the ID na¬ 
ture of the universal eigenchannel structure. For diffu¬ 
sive samples without edge reflections, Fi{x) has an ex¬ 
plicit form of 'KxiyL — x)l{2Ltj which is symmetric with 
respect to the sample center (Fig. red curve in up¬ 
per panel). This quadratic form of Fi{x) can be ob¬ 
tained by solving the (normal) diffusion equation with 


the Dirichlet boundary conditioiP. In the presence of 
edge reflections, the diffusion equation is the same but 
the boundary conditions are of the mixed type. Such 
diffusion equation can be easily solved, giving Fi{x) = 
tt{x + — X Zb)/{2{L + z^ -t - Zb)t). Upon substi¬ 

tuting this expression into Eq. ( [m| we find (Fig. blue 
curve in upper panel) 


Wr=i{x) = l + 


Tr{x + zl){L - X + Zb) 
2{L z'^-\- Zb)i 


(55) 


When Zb, vanish Eq. ( [^ reduces to the result obtained 
in Ref. H In Appendix |C| we further discuss a relation 
between this eigenchannel profile and the field theory (|^ . 
The position of the center of the profile (551 is pre¬ 


cisely the same as what is given by Eq. (46). That is, for 


this eigenchannel (T = I) the resonator center is locked 
to the profile center. So, if the reflection of the right 
(left) edge is stronger than that at the left (right), i.e., 
Zb > {zb < z'ij) the profile and resonator centers move 
to the right (left) off the sample center (cf. Fig. [^. As 
the asymmetry parameter continues increasing, eventu¬ 
ally the resonator center approaches the sample edge, and 
the perfectly transmitting eigenchannel disappears, sig¬ 
naling the DTE transition. We emphasize that although 
the resonator is due to the interference of multiply scat¬ 
tered random fields, it exists in the diffusive sample, far 
from the localization regime. 


VI. CONCLUSION AND DISCUSSION 

We have shown analytically and confirmed numerically 
that asymmetry in the edge reflections has a significant 
impact on the wave propagation through disordered me¬ 
dia. In particular, it can trigger a peculiar DTE double¬ 
transition phenomenon in quasi ID diffusive samples. 
When the asymmetry in the reflections of the two sample 
edges is weak so that the inequality |C ~ Cl < 1 is satis¬ 
fied (i.e, the difference between the two edge resistances 
is smaller than the bulk resistance), high transmission 
eigenchannels are present, and the DTE exhibits a sin¬ 
gularity p{T —>■ I) ~ (1 — T)“5. In the opposite case 
of strong asymmetry (i.e, the difference between the two 
edge resistances is larger than the bulk resistance), high 
transmission eigenchannels are closed. When the asym¬ 
metry parameter |C ~ C^l = 1 (i-6) the difference between 
the two edge resistances equals to the bulk resistance), 
the system is in a critical phase with the DTE transition 
exhibiting critical statistics, i.e., p{T —>■ I) ^ (I — T)“3. 
These phenomena are universal in the sense that they are 
governed by a single parameter | C ~ C^ I ^md at each phase 
regime, the asymptotic behavior of p{T —>■ I) does not 
depend on the details of system’s structure such as the 
specific value of edge reflection, disorder configuration, 
sample width, etc.. 

Surprisingly, notwithstanding its occurrence in diffu¬ 
sive samples, far from the localization regime, the DTE 
transition has a close connection to the resonator model 
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of strong localization in ID. More precisely, the DTE can 
be mapped onto the statistics of the resonator center so 
that the factor /(= p/po) corresponds to the spatial den¬ 
sity of resonators. We show that perfectly transmitting 
eigenchannels can be modeled by ID resonators, and the 
disappearance of such an eigenchannel as the asymmetry 
in the edge reflections increases mimics the evolution of 
a transmittance resonance when the corresponding effec¬ 
tive cavity shifts to an edge of the sample. 

Our findings indicate a novel coherent phenomenon of 
diffusive waves. They exist in very general wave systems, 
e.g., quantum matter and classical elastic waves, wher¬ 
ever edge reflection is strong. A prominent system is the 
normal-metal - superconductor junctiorP^, where edge 
reflection is created by the tunnel barrier at the normal- 
metal - superconductor interface. Another system is the 
dissimilar solicP^ where the Kapitza resistance on the 
thermal phonon transport are not negligible. Our find¬ 
ings may enable the control of transmission eigenchannels 
and eigenvalues via edge reflection, which is well within 
experimental reach. The indication of the existence of 
resonators in diffusive samples may find practical appli¬ 
cations such as low-threshold lasing. 

An important subject of further studies in this direc¬ 
tion is to better understand the formation of resonators 
in quasi ID diffusive random media. A related issue is 
the connection between the resonator and the eigenchan¬ 
nel structure for arbitrary transmission eigenvalue. This 
would help to uncover the physical meaning of the lin¬ 


ear connection (42) between the parameter <f> and the 
resonator center x. The subject of particular interests 
is the interplay between edge reflection and gain (or ab¬ 
sorption) . 

Our investigations of the DTE transition have been 
restricted to quasi ID diffusive samples. We should em¬ 
phasize that the microscopic formalism, i.e., the super- 
symmetry field theory of DTE presented in Sec. m can 
be directly applied to localized samples. The additional 
technical difficulty in calculating Eq. ([^ is that one needs 
to take into account all the saddle point configurations 
in which the supersymmetry is broken. How would the 
DTE transition be affected then? We leave this challeng¬ 
ing problem for future studies. 
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Appendix A: Derivations of Eqs. (29| and (30) 


We rewrite Eq. (27) as 


2CI + ACl/2 


sinh^ '0+ (a cosh + b) + sinh^ 'ip_ (a cosh + h) 

(a cosh 0+ -I- b)(a cosh 0_ -|- b) 

a(sinh^ 0_i_ cosh'0_ -|- sinh^ 0_ cosh0_|_) -|- 6(sinh^ 0_|_ -|- sinh^ ip_) 
a? cosh'0_|_ coshi/)- -|- a&(cosh'0+ -|- cosh^-) -|- b'^ 


and Eq. (28) as 




sinh^ 0+ (a cosh tp- + b) — sinh^ 0_ (a cosh 0+ -|- b) 

{a cosh 0+ -I- b) {a cosh ip- + b) 

a(sinh^ 0+ cosh0_ — sinh^ ip- cosh'0+) -|- 6(sinh^ 0+ — sinh^ ip-) 
cosh 0_|_ cosh ip- + a&(cosh ip+ + cosh ip-) +b‘^ 


(Al) 


(A2) 


To proceed further we use the following identities: 


cosh 0+ cosh 0_ 
cosh 0+ -I- cosh 0_ 

cosh 0+ — cosh 0_ 


(cosh 20 -|- cosh A0)/2, 
A0 

2 cosh 0 cosh , 

- . Alp 
2 sinh Ip sinh , 


with Ip = (0_|_ -I- 0_)/2 = C^- (p, Alp = 0_,_ - 0_ = 
AC^ — 2iTT, to obtain 

sinh^ 0+ cosh ip- + sinh^ ip- cosh 0+ 

= (cosh 20-|-cosh A0 — 2) cosh 0 cosh(A3) 


sinh^ 0+ cosh ip- — sinh^ ip- cosh 0+ 

= (cosh20-|-cosh A0-I-2) sinh0sinh —(A4) 
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sinh ?/;+ + sinh ip- = cosh 2-!/i cosh A'i/; — 1, (AS) Substituting them into Eqs. (Al) and (A2) gives 


and 


sinh^ -0+ — sinh^ '0_ = sinh 2ip sinh /S.ip. (A6) 


a(cosh 20 + cosh A0 — 2) cosh 0 cosh + 6(cosh 20 cosh A0 — 1) 

2C, + ACj/2 = ---- i— - 

(cosh 20 + cosh A0)/2 + 2ab cosh 0 cosh ——h b'^ 


(A7) 


2C.AC'^ = ^ 


a(cosh 20 + cosh A0 + 2) sinh 0 sinh + 6 sinh 20 sinh A0 


(cosh 20 + cosh A0)/2 + 2ab cosh 0 cosh + 5^ 


(AS) 


which are Eqs. (29) and (30) upon the substitutions of AC^, 0, and A0. 


Appendix B: Values of a and /3 


Taking into account 6 — a = 1 and introducing / = tt/, we rewrite Eqs. (29) and (30) as 


and 


^ 2 , 


2 ( 0 ^ - / ) = 


4(0 + 0)/ = 


—a(cosh20 + cos 2/ — 2) cosh 0 cos / + (a + 1) (cosh 20 cos 2/ — 1) 
a2(cosh20 + cos 2/)/2 — 2a{a + 1) cosh0 cos / + (a + 1)^ 

—a(cosh 20 + cos 2/ + 2) sinh 0 sin / + (a + 1) sinh 20 sin 2/ 


a^(cosh20 + cos2/)/2 — 2a(a + 1) cosh0cos/ + (a + 1)^ 


(Bl) 


(B2) 


As shown in Sec. HIE 1[ the numerical solution to 


Eqs. (29) and (30) yields the general form, 


/ = fci0“, 0 = k2(, 


0 —>■ 0 , 


(B3) 


where fci 2 are coefficients depending on a. The numerical 
solution further yields a (3 1/3. 


We now show that Eq. (B3) and the relation above for 


a and (3 lead to a stronger result. Let us expand Eq. (Bl) 


around 0 = 0. B ecau se of Eq. (B3), the first few order 
terms of 0 in Eq. (Bl) are 00, / , 0^/ and / . At this 


stage we cannot determine for these terms which one is 
smaller, but all the other terms in the expansion must 
be of higher order because o f a » j3 ~ 1/3. We perform 
the same expansion for Eq. (B2), for which the first few 


_ _2 — _ 

order terms are 0/, / 0 and /0^. As a result, 

400 - c(7^ - 6/^0^ + 0'‘) = 0, 
40/ + 4c(7^0 - 70^) = 0, 


(B4) 

(B5) 


where the coeffi cient c = a + 2/3. 

Because Eq. (B5) is satisfied for arbitrarily small 0, 


we immediately find a = j3 = 1/3. However, to justify 


Eq. (38) we must further prove that are rea l. To this 
end we su bstit ute a = /3 = 1/3 as well as Eq. (B3) into 
Eqs. (B4) and (B5), obtaining 


4^2 — c{k\ — &k\k\ + fc^) = 0, 
4fci + Ac{k\k2 — kik2) = 0. 


(B6) 

(B7) 


Upon eliminating the coefficient c we reduce them to kf — 
2 kfk 2 — 3^2 = 0 which gives 


kl - 3fc| = 0. 


(B8) 


Substituting it into Eq. (BIO) we obtain ki = 
3^/^(2c)“^/^, fc 2 = —(2c)“^/^. As a result, 

7 = 3i/2(2a + 4/3)-i/30i/7 (B9) 

0 =-(2a + 4/3)-i/30i/^ (BIO) 

Summarizing, combined with numerical analysis we have 
shown that Eqs. (B9) and (BIO) are solutions to Eqs. (29) 
and dSOl). 
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Appendix C: Relation between field theory and 
universal eigenchannel structure 


In this Appendix we show that in the functional in¬ 
tegral formalism, fluctuations around the saddle point 
carry information on the universal eigenchannel struc¬ 
ture. To illustrate this we focus on the perfectly transmit¬ 
ting eigenchannel, i.e., the transmission eigenvalue 7^=1 
or equivalently (j) = Q. Because of this the functional in¬ 
tegral is reduced to 


r(2ziQd^Q-[Q,K\)U=L=0 


' (2z'^Qd^Q+[Q,K])U=o=^ 


Zl[g](-)e"t/o dxst-cid^Q) 


Here the pre-exponential factor (•) depends on specific 
observable considered. Its details are unimportant for 
present discussions of the general structure of field theory. 
Comparing with Eq. Q, most importantly, the boun dary 
constraint of the right sample edge is replaced in Eq. (Cl) 
by {2zbQd^Q - [Q,I^])\x=l = 0 (noting that = 

A) while that on the left remains the same. 

To proceed further, we introduce the rational 
parametrization, 


Q = {l + iW)K{l + iW) 


-1 


(C2) 


where the supermatrix W = anticommutes with 

A, implying the ar-sector index A ^ A'. Recall that a, a' 
are fb-sector indexes. Then, we expand Q in terms of W 
and substitute the expansion into both the action and 
the boundary constraints (as well as the pre-exponential 
factor). Eor diffusive samples {L we may keep the 

expansi ons up to the quadratic order. As a result, we 
reduce ( |C1[ ) to 

/ D[W\{-)e-i^o dxstr(d^W) ^ ((^3) 

7(z'9,-i)iyu=o=o 


where we have used the fact that the Jacobian for the 
transformation: Q —>■ IT is unitjEH. 


The propagator of the effective field theory (C31 is none 


but the well-known diffuson in diagrammatic techniques 
(recall that we do not consider time-reversal symmetry 
in the present work). More precisely, this propagator is 
given bjim 


= ^yo{x,x'). (C4) 




1 


I.e., 


/ /^ 1 {z'b + x^){L + Zb — X:^) 

yo[x,x} = — - . , . , ,, -, (C6) 


Dn 


■ Zb ■ 


where = min (max) {a;, a;'}. 

Physically, yo{x,x') is the intensity (energy density) 
profile when a unit flux is injected at x'. The expres¬ 
sion (C6| entails the canonical meaning of Zb^z'b - the 
so-called extrapolation lengtlP^. As shown in Fig. [T^ 
the profile linearly falls down from the source point and 
vanishes outside the sample and at a distance of z'f^ (zb) 
to the left (right) end. In the limiting case of vanish¬ 
ing Zb, z'l^, as shown in Ref. SI 3 ^ 0 ( 2 ^! x' = x) is essentially 
El (x) in Eq. (54) up to an irrelevant overall factor and 


thereby determines the eigenchannel profile WV=i(a;)- If 
we assume that the expression (54) is valid also for non¬ 


vanishing Zb,z'f^, then Eq. (55) is received. 



FIG. 12: The energy density profile, yo{x,x'), generated by 
injection of a unit flux at point x' inside the sample. 


Appendix D: Discussions on parametrization (42) 


In this Appendix we show that provided x is uniformly 
distributed in the range of [0,L], Eq. (42) is the unique 
parametrization leading to the bimodal distribution. For 
the convenience below we introduce the variable 77 = xjL. 
We have 


dr] a' 

^ ~ Tv/TAAT’ 


(Dl) 


where the coefficient a' is fixed by the normalization con¬ 
dition: 



ryrATr 


dT = l. 


(D2) 


Here (-jo sta nds for the average with respect to the Gaus¬ 
sian weight ( |C3| ) (subject to the boundary constraints), 
and V is the average density of states. The localization 
length ^ = 2'kvDq with Dq being the Boltzmann diffusion 
constant. Note that propagators with other index com¬ 
binations vanish. Importantly, yQ(x,x') is the solution 
to the normal diffusion equation: 

- DQdlyQ{x,x) = 5{x-x), (C5) 

{z'bdx - l)yoix,x')\x=o = 0 , 

( Zbdx + l ) yo { x , x ')\ x=L = 0 , 


Since T{r]) is symmetric with respect to the middle of 
the samp le, 77 is a double-valued function of T. Solving 
Eq. (|DI[) gives 


_ / 1 — yi —^ 1 + yi ~ Ttyp 

\l + ^qrzrri-^l-Ttyp 


(D3) 


In combination with Eq. (D2) it gives 



In 


1 + ~ 'Tt , 


typ 




typ 


± In 


1-VT^ 
1 + VT^ 


(D4) 
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where the + (—) sign corresponds to rj < 1/2 (> 1/2). 
From this expression we find 


r 


b' 



1 + t/1 — Tty 


(D5) 


Substituting Ttyp = [cf. Eq. (51)] into the ex¬ 

pression of b' we find b' « 2L/£. Equation (D5) then is 
identical to Eq. (42). 
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